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Abstract 

We propose a realistic model with Majorana neutrinos in the framework of uni¬ 
fying the three generations of fermions by point interactions in an extra dimension. 
This model can simultaneously explain the origin of fermion generations, fermion 
masses and mixing, and the smallness of the masses of Majorana neutrinos. We show 
that there are two mechanisms working together to suppress the neutrino masses 
significantly, so we do not have to introduce a very large extra-dimension cut-off 
scale. One is the type-I seesaw mechanism and the other is the overlap integration 
of localized lepton wave functions. A singlet scalar with an exponential-like VEV 
plays a central role in these two mechanisms. For consistency in this model we intro¬ 
duce a U{1)' gauge symmetry, which will be broken by the singlet scalar. Parameters 
of our model can fit the masses and flavor mixing data well. These parameters can 
also predict all CP violating phases including the Majorana ones and accidentally 
rescue the proton from decay. 


1 Email: zhh98@mail.sysu.edu.cn 


1 



Contents 


1 Introduction 

2 The model 

3 The lepton sector 

3.1 1/(1)' symmetry and type-I seesaw. 

3.2 Numerical results and discussion. 

4 Summary 

A The general setup of the framework 
B Quark masses hierarchy and flavor mixings 
C Why an explicit Majorana mass term does not work 


a 

s 


17 

18 

25 

28 


1 Introduction 

The recent discovery of the Higgs boson is a great success for the Standard Model 
(SM) of particle physics. In the SM, the masses of weak gauge bosons and fermions 
are generated by the Higgs mechanism, which predicts the existence of a CP-even 
scalar particle, and finally this only scalar boson was discovered at the Large Hadron 
Collider (LHC) in 2012 [3, 0 • 

However, many people believe that the SM should not be the finale of particle 
physics. One of the reasons is that it cannot explain the large hierarchy of fermion 
masses. In the SM, all fermion masses, mixing angles and CP phases are free pa¬ 
rameters. If one looks at the mass spectrum of fermions, one will find a significant 
hierarchy between different generations. The hierarchy between quark sector and 
lepton sector is even worse. 

In the original version of the SM, the neutrino masses are assumed to be zero. 
However, to explain the oscillation phenomena observed in experiments, the neutri¬ 
nos have to be massive. Similar to the way used in the SM to give fermions masses, 
it can make neutrinos massive by introducing right-handed neutrinos which couple 
to the Higgs field through Yukawa terms. But this way is quite unnatural due to the 
large hierarchy. A cosmological observation from Planck set a 0.23 eV upper bound 
for the sum of the three generations of neutrinos [ 3 J. It leads to about 11 order of 
magnitude hierarchy between the Yukawa coupling of top quark and the neutrinos. 
This unnaturalness indicates us a strong motivation to go beyond the SM. 
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There are three types of seesaw mechanisms to explain the smallness of neutrino 
masses. The type-I seesaw introduces right-handed neutrinos coupled with the left- 
handed leptons through Yukawa interactions, and then the Majorana masses of the 
left-handed neutrinos will be generated by a higher dimensional operator and be 


suppressed by the heavy Majorana masses of the right-handed ones [4H6j]. The type- 
11 seesaw introduces triplet scalars coupled with the left-handed lepton doublets, 
and the vacuum expectation value (VEV) of the scalar will be suppressed by its 
large quadratic masses [3-0] . The type-111 seesaw is similar to the type-I, but it 


introduces heavy triplet leptons 10]. All these mechanisms usually need a high 


seesaw energy scale, for example the Grand Unification Theory (GUT) scale, to 
suppress the induced Majorana masses of the left-handed neutrinos. 

Besides the seesaw mechanisms, an alternative way to explain the masses hi¬ 
erarchy naturally is to enlarge the spacetime dimension. One interesting case is 
the thick wall model 11 , in which fermions have Gaussian wave functions of the 


5th dimension coordinate and their locations are determined by their 5-dimensional 
(5D) masses. When two fermion wave functions are separated slightly, their over¬ 
lap integration with the Higgs VEV profile will be suppressed exponentially, then 
a large hierarchy structure between fermions can be naturally obtained. Another 
fascinating case is the Randall-Sundrum model 12, Il3| . in which right-handed neu¬ 
trinos localize near a hidden brane, while the other fermions and the Higgs held are 
confined on a visible brane. Thus the right-handed neutrinos interact with the other 
fields weakly, and they only have tiny masses. 

Recently, a new extra-dimension model Ml, [l5j was proposed to unify the 3 
fermion generations. The model introduces 5D fermion fields living in an extra- 
dimensional interval or circle with several point interactions (he. 0-thickness branes). 
For each 5D fermion, there are three independent modes between branes. They be¬ 
have as three generations, and the hierarchy between generations is achieved by 
coupling the 5D fermion held to a scalar held which has an exponentially increasing 
extra-dimensional coordinate-dependent VEV. This specihc VEV can be generated 
by imposing Robin’s boundary conditions on the scalar at two boundaries of the 


5th dimension (see more details on the phase structures in 16]). In addition, a 


twisted boundary condition is imposed on the Higgs doublet to create CP violating 
phases for both quark and lepton sectors jT3] • In Ref. [10], a 5D singlet neutrino 
held (which has a right-handed chiral neutrino 0-mode) is introduced to construct 
Dirac masses terms for neutrinos, and the smallness of neutrino masses are obtained 
from a proper arrangement of the point interaction positions. A stringent constraint 
on the model with a set of htted parameters is to suppress the proton decay rates. 
By a rough analysis with some baryon number violating dimension-8 operators, the 
cut-off A ~ L" 1 is estimated to be as large as the GUT scale (10 15 GeV). 
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In this paper, we discuss a possibility to extend the model of Ref. [l5| to a Majo- 
rana neutrino case and to avoid the large cut-off scale. To implement this, we need 
a Majorana mass term of the singlet neutrinos. A naive trial is to write down an 
explicit Majorana mass term for the singlet neutrino fields and their charge conju¬ 
gation. However, it fails since the equations of motion for the singlet neutrinos no 
longer respects the so called quantum mechanical supersymmetry (QMSUSY) which 


is important for acquiring chiral zero modes ll|, |l8|. The existence of a Majorana 


mass term implies that the lepton number is no longer a conserved quantity, and 
thus a dimension-7 effective operator La 2 H* a 2 L c may appear in the Lagrangian 
in principle. Here L(x,y) is the 5D lepton doublet held, H(x,y) is the 5D Higgs 
doublet, and the power counting is achieved in 5D spacetime. But this effective 
operator can induce large Majorana masses of the left-handed neutrinos after the 
electroweak symmetry breaking. To avoid large neutrino masses which violates the 
experimental bounds, it requires either a high cut-off scale or a very small coupling 
constant for this term. 

To overcome this problem and to forbid the harmful explicit Majorana masses 
terms at the same time, we introduce a new t/(l)' gauge symmetry. If we let the 
singlet neutrino held Nr and the combination H^a 2 L c be U(l)' charged, none of 
these two annoying terms, N^Nr and La 2 H*H^cr 2 L c , can survive, since they are 
doubly U{1)' charged. In other words, we increase the symmetry of the model to 
prohibit the unwanted operators like La 2 H*a 2 L c . 

We need a natural way to realize the experimentally acceptable Majorana masses 
for neutrinos. Consider another dimension-7 operator ($*) 2 N^Nr + h.c., which is 
available if the U (1)' charge of the singlet scalar $ is assigned to be the same as that 
of the singlet neutrino Nr. Obviously, This term can contribute to a Majorana mass 
for the right-handed neutrino 0-mode when the scalar $ obtains a non-zero VEV 
and break the f/(l)' gauge symmetry. The 5D scalar <h(x, y) is initially introduced to 
realize the hierarchy of the three generations of quarks and leptons, and it is imposed 
on the Robin’s boundary condition to get a VEV, (&(y)), as an exponential-like 
function of the extra dimensional coordinate y. This VEV has the effect of killing 
two birds with one stone. If the O-thickness branes’ positions of singlet neutrino are 
chosen appropriately, that is, if the third generation singlet neutrino wave function 
has a big overlap with the large value side of ($(?/)), it can obtain a mass which 
is much larger than the Dirac masses. This large mass turns on the type-I seesaw 
mechanism to lower the neutrino masses further. 

To make our model self-consistent, we set the U{1)' charge of each held agreeing 


with the anomaly free conditions 19]. We also consider the constraint from the pro¬ 


ton decay. By some simple analysis with the dimension-8 baryon number violating 
operators, we see that for our best-fit parameters the proton will not decay. So it is 
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not necessary to let the cut-off energy be the GUT scale in this model. 

An outline of the paper is as follows. In section [21 we will summarize some key 
elements of the model and building a realistic model in the framework. We also 
discuss what are the problems of introducing an explicit Majorana mass term. In 
section [3], we discuss how to generalize the model to include Majorana neutrinos and 
how the seesaw mechanism works with a few TeV extra-dimension energy scale. We 
will also fit the data of leptons and do some discussion. Section m is a summary. In 
Appendix lAllBIICl we supply some mathematical details of the discussion in section 

El 


2 The model 


To begin with, let us summarize some general setups of this model (with some 
mathematical details reviewed in Appendix [AJ) HQ: 


• The spacetime is extended by a finite size of space-like extra-dimension, i.e. 
an interval or a circle. The mode expansion is made as usual and the lowest 
modes, i.e., the zero modes, are regarded as the SM particles. The mass gap 
between the 1st. K.K. modes and the zero modes is roughly the inverse of the 
5th dimensional size. In many extra dimension models, the mass scale is at 
least around the energy scale of the LHC experiment. 


In the free field limit, there is a quantum mechanical supersymmetry (QM- 
SUSY) between the left-handed and right-handed components of 5D fermion 
[ill . H. 221. This symmetry ensures that the left-handed and right-handed 
modes at the same level have equal masses. Thus, their 4D parts can be sep¬ 
arated from the 5th-dimension-coordinate-dependent parts, and can form a 
Dirac fermion satisfying the 4D Dirac equation. In particular, for the zero 
mode, the symmetry together with the Dirichlet boundary conditions implies 
that one of the chiral spinors should vanish and the other one is massless. This 
is the method of generating chiral zero modes. 


An important ingredient for unifying generations is the point interaction [14, 


15j, which can be regarded as a Delta-function-like interaction. This specific 
interaction is located at a point in the 5th dimension and results in the Dirichlet 
boundary condition for the 5D fermion. If we introduce two interacting points, 
then they will separate the interval at extra dimension into three pieces. The 
modes living in different pieces are independent from each other although they 
come from the same 5D fermion held. These different modes can be regarded 
as different generations. 
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To achieve the hierarchy among generations, a singlet scalar held $(#,?/) is 
introduced to couple with 5D fermions. A Robin’s boundary condition on the 
5D scalar will force its VEV ($(?/)) to be ^-dependent as 


<*(y)> = 


cn (v 2k(y~yo)’ k ) 


(i) 


where the function cn(x, k) is the Jacobi elliptic function of x with index k , 
and k,y,v are dehned as 


,2 


fl 2 +U 2 

/i 2 =¥(l + 

M 2 


- 2 = w(\A + tf 


-i I 4A|Q| \ 
1 ~r~ Af 4 > 

4A|Q| 


( 2 ) 


A stud y o f this 
An 
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with Q,yo being constants of integration determined by L±. 
singlet scalar with Robin’s boundary condition can be found in Ref. 
important result in their study is that <E>(x, y) can couple with gauge fields 
corresponding to some group, such as a U(l)' group. This symmetry will break 
if L < L c = pb_tanh _1 ( ^ — > 1S|- Usually we use the condition 
M 2 < t- 2—,T max = max(L. ,L_), which is sufficient but not necessary. 

■L* max x 7 

When we proceed to construct a realistic model comparable with experiments, 

The requirements are briefly listed 


some special settings are also needed [14|, lla 117 
as follows: 


1. The 5th dimension need to be a circle (S' 1 ). This is a part of the requirements 
from the flavour mixing behavior of the SM. And it is also consistent with the 
twisted boundary condition setting of the Higgs doublet. 


2. We need to specify the 5D matter fields with appropriate boundary condi¬ 
tions. In the quark sector, we should introduce an electroweak SU (2) doublet 
quark Q(x,y ) = ( U^x^y ) D L (x,y)) T , and two singlets quarks Un(x,y) and 
Dn(x,y). For the doublet Q, we use a Dirichlet boundary condition PrQ = 0 
at y = L = 0, L^\l^ so that its zero modes are left-handed, while for 
the singlets Ur and Dr, we use Dirichlet boundary conditions PlUr = 0 at 
y = Lq U \ L[ a \ and PlDr — 0 at y — L ( f f\ L ( f \ so that their zero 
modes are right-handed. Note that in general L^ q> are different from and 
iJf 1 . This is necessary for flavor mixing structure. For the lepton sector, the 
situation is similar to the quark case. We just replace the quark doublet by 
a lepton doublet and the up and down type quark singlet by neutrino and 
charged lepton singlet. 
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3. We need a Higgs doublet H(x, y) to couple with fermion fields through Yukawa 
couplings. Of course it should acquire non-zero VEV (H) to break the elec- 
troweak symmetry. A special treatment is to impose a twisted boundary con¬ 
dition on H(x, y ) as H(y + L) = e ld H(y) [l7]. This twisted boundary condition 
will make the VEV (H) get y dependent phase as ( H(y )) = -^=e^ y , then its 
overlap integration with fermions’ wave functions will produce CP phases for 
CKM or PMNS matrices. 

As an example, the detailed treatment of the quark sector are presented in Appendix 
iBl We also fit the parameters of quark sector independently and list them in Table. [4] 
The fitting will fix the M parameter from the singlet scalar $ and the 6 from the 
Higgs H , and they will be regarded as input data for the lepton case. 

Before going to the next section to discuss our treatment of the lepton sector. It 
will be helpful to ask what’s wrong if we just write down an explicit Majorana mass 
term? We will discuss this briefly as follows, and supply more details in Appendix 
□ 

One problem of this naive trial is that Majorana mass term will modify the 
equation of the motion for the 5D fermion. This modification breaks the QMSUSY 
between the left-handed and right-handed components in the E.O.M. As we men¬ 
tioned previously, generating chiral zero modes rely on this symmetry. 

Another problem with this naive trial is that since we are going to break the 
lepton number conservation explicitly, then in principle we should also include an 
operator as La 2 H* a 2 L c which has the same dimension with the terms we used 
to generate the Dirac masses for leptons. After the Higgs acquires a non-zero VEV, 
this operator will generate Majorana masses for the left-handed neutrino zero modes. 
Then a fine-tuning is needed when we diagonalize the neutrino mass matrix to obtain 
sub-eV masses. 


3 The lepton sector 

3.1 U(l)' symmetry and type-I seesaw 

For the lepton sector, we introduce an 577(2) doublet L = (N L (x,y), E L (x,y)) T , 
and singlets N R (x,y), E R (x,y). When we consider the structure of our model, the 
lepton number is not necessary to be preserved. The most famous model which 
violate lepton number is the type-I seesaw [7j. In type-I seesaw a Majorana mass 
term for the right-handed neutrino is introduced. If the Majorana mass M R is 
extremely large comparing to the Dirac mass , then after diagonalize the mass 
matrix, a mass for the three lightest neutrinos taking the form —m^M^ 1 m^ T 
will be suppressed significantly. But as we discussed in the section [2] an explicit 
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Majorana mass term is not allowed to exist. We will assign a U( 1)' charge to Nr to 
forbid such a troublesome term to keep the chiral 0-mode, and then use the VEV of 
the scalar <f> to create the Majorana masses for the right-handed neutrino 0-mode. 

As we have mentioned in section [2j the La 2 H* a 2 L c operator will bring us 
a problem of fine-tuning. To solve this problem it will be forbidden by the U( 1)' 
symmetry if we let Lia 2 H* to be charged. All these indicate that we would better 
add the U{1)' symmetry into the model. Then to justify the model, we should put 
some constraints to the undetermined U(l)' charges. 

The gauge group in our model is now SU(3)c x SU(2) L x U( l)y x U( 1)'. Let’s 
denote the representation of all left-handed zero modes in the form (N Cji , N W}i , Y j, Q 7 ), 
where N Cji and N Wii denote the dimensions of SU (3) and SU (2) representation(conjugated 
representation with a bar) of f-tli field, while Y t and Q\ denote the U( 1) hyper-charge 
and U{1)' charge of i-th held. N Cji , N W}i ,Yi for each type of held are just the same 
as iu the Standard Model. Q \s for each type of held are unknown variables and will 
be determined later. We list the representations for fermions in table |T] Now the 


Table 1: Gauge group representations for fermions 


Fields 

q 

U R 

d C R 

Reps. 

(3, 2,1/6, Q' g ) 

(3,1, —2/3, Q' u ) 

(3,1,1/3 ,Q' d ) 

Fields 

l 

V R 

e R 

Reps. 

(1,2,-1/2, QO 

(1,1,0, %) 

(1, 1, 1, Q'e) 


covariant derivatives for each held are 

D^ ] = d N - ig s G l N t l - igW^T a - i\g'B N - iQ' q g c C N 
D% ) = d N - ig s G l N t l - i\g'B N + iQ' u g c C N 
= djsr — ig s G l N t l + i^g'Bjy + iQ' D g c C n 
= d N - igW^T a + i\g'B N - iQ\g c C N 
' = d n + iQ' n g c C]\r 
D^ ] = <9v + ig’B N + iQ' e g c CN 

= d N - igWyjT a - i\g'B N - iQ' h g c C N 
D^ n ^ = On — iQ'^cCN 


Where the CV is the gauge held corresponding to U( l) 7 and g c is the gauge coupling. 


There are 6 constraints of Q\ come from the consideration of anomaly free 19|. They 












are 


2 Q'q + Q'u + Q'd — o 

3g; + Q\ = 0 

e(|) 2 g; + 3[-(|) 2 g; + (±) 2 g;] + 2 (-|) 2 g; + g; = o 
eg? + 3[g; 3 + q '!] + 2 g; 3 + g; 3 + g; 3 = o 
6 • \ q ' 2 + 3[-§g; 2 + \ q ' 2 } + 2(-i)g; 2 + g' e 2 = o 
eg; + 3[g; + g;] + c iq \ + g; + g; = o 


( 4 ) 


It seems that we have 6 equations for 6 variable, but actually only 4 of them are 
independent. We rewrite g's in terms of Q\ and Q' e as follows 


{ Q q = -Wi 
Q'u = -I 
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Q'l - Q'e 
Q' d = iQl + Q'e 
Q'u = -2 Q'l - Q'e 


(5) 


Then when we choose a set ( Q[, Q ' e ), all the other variables are determined. For our 
purpose, we will impose more theoretical constraints on g's. One is that we need 
Yukawa terms as 


$Q(ia 2 H*)U R , $*QHD Rl 
$L(ia 2 H*)N R , <&*LHE r . 

to be gauge invariant. Assign a U(l)' charge Q' h to H and Q', to $, and use 
finally we find the only constraint is 


( 6 ) 


Q'i + Q'e ~ Q'h + Q'<j> — o 


(7) 


Another important constraint is to let Q^ = Q' n so that ^ 2 N^N R is gauge invariant, 
or let g; = —Q' n so that ^* 2 N^N R is gauge invariant. Then we replace Q ^ by 
ig; in (ED and use (ED, we obtain Q' h = —Q\ for ^ 2 N^N R or Q' h = 3g^ + 2 Q' E for 
Q* 2 NftN R . Remember that we want L(ia 2 H*) to be 1/(1)' charged and it requires 
that Q' h 4 —Q'i, so only Q' h = 3g( + 2 Q' e corresponding to ^* 2 N^N R is allowed. Of 
course, we should have Q' n ^ 0 to kill the explicit. Majorana mass term for singlet 
neutrino and this requires that Q' e ^ —2Q[ . The other constraints may come from 
experimental considerations but that is beyond the scope of this article. 

There are still many possible choices of g's and we only list three interesting 


candidates which are similar to 23, [29 


1. U R \ Q\ = Q' q = 0, Q' u = 1, Q' d = -1, Q' e = -1, Q'n = 1, Q' h = -2, g' 0 

2. u B — L \ g; = u = Q' d = 4,g; = -i,g; = g; = i,g; = -i,g; 


-i. 

= -i. 
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3 . U x . Q' q = \,Q' U = \ )Q ' d = ~l,Qi = -l,Q' n = l,Q'e = \,Q'h = -I’Q* = “ 1 - 

The mass term of zero-mode leptons will be generated by 


Cyuk = ~ J dyiy^$(y)L(ia 2 H*)N R + y^$*(y)LHE R + h.c.} 
-I [ dy[y^^ 2 W R N R + h.c.} 


( 8 ) 


where y^ n \yN) and ?/ m ) are couplings with dimension —2. After the U( 1)' and 
SU( 2) x U( 1) breaking, two terms in the first line generate Dirac mass matrices for 
charged leptons and neutrinos and the term in the second line generate a Majorana 
mass matrix for right-handed neutrinos. 

Imposing Dirichlet boundary conditions on fermion fields, twisted boundary con¬ 
dition on Higgs doublet and Robin boundary condition on $, we can expand fields 
in modes and finally obtain their profiles: 

L = E / ( \ (0)/ J + ( KK modes), 

3 

( 0 )/ 


Er = 52 f e ( 3 )(y)ej$(x) + (KK modes), 


2—1 

3 


Nr = 


E fMy)^R( x ) + ( KK modes), n r = CN r 

Z -* u iR 

2 — 1 

/,«(») = e(B-h- 1 )9(h ,) -!/). 

L iL 

/>(») = N§e- M ^- L ^e(y - iftWd'’ - y), 

iR 

/_«(») = - d-iMd 0 - y) 


(9) 


where N - 1 ^, are normalization constants. Substituting these profiles into 

OSD, we get the Dirac mass matrices and Majorana mass matrix: 


= Jdyy^^.my))f l m(y)f e m(y)e'L, 
m( D,ij = I d yy {n) ^r(^(y))f l w(y)f u m(y)e~^, 


iSy 


r(y) 


iSy 


( 10 ) 


Mr, a = y (m) J,m %(*('! i)) 2 f„m(y)f„m(y) 


in 


Obviously, Mr is a diagonal matrix since the integration only involves the profile of 
Nr. Now we write the chiral zero modes in Weyl basis: 


..(0) v .. jo) . . 

U iL U iL,a, Oi e iL,a, 


(0) f ,a 

V iR U iR 


( 11 ) 
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where a, a are indices of Weyl spinors. Then for neutrinos we can represent the mass 
term as 


r 0) — 

^mass 


2 \ V iL,a U iR,d 



+ h.c. 


( 12 ) 


Following Xing’s parametrization and discussion 24], we introduce a 6 x 6 unitary 
matrix U to transform the mass eigenstates to flavor states. U can be decomposed 
into 


(1 0 \ (A R\ (Vo 0 \ = (AV 0 R\ 
VO Uo) \S BJ V 0 i) \U 0 SVo UqB J 


(13) 


where Vo and Uq are 3x3 unitary matrices and A, B, R, S are 3 matrices under the 
unitary conditions: 

AA f + Rtf = BB t + S& = 1, 

AS* + RB t = ARi + S^B = 0, (14) 

A f A + &S = B^B + RtR = 1 


We can use U to diagonalize the mass matrix in (TT2|) : 




(( m Ap 


)T 


( n) \ 

m D,ij ] W 

M R,ij ) 


(% jo\ 
V o m n J 


(15) 


where M v and Mjy are diagonal matrices: M u = Diag{mi,m 2 ,m 3 } are very small 
while M n = Diag{Mi, M 2 , M 3 } should be very large. Finally we can hnd approxi¬ 
mately 


M v ~ -Vl(m%-Ml x {m 


( n ) ) 
D,ij) 


> 0 * 


(16) 


The minus sign can be absorbed into charged lepton basis. Remember that at the 
beginning of this section, we use the {7(1)' symmetry to kill the La 2 H*H^a 2 L c 
dimension-7 operator. When the U{1)' symmetry breaks spontanuously, this term 
comes back by connecting two Yukawa interaction with an internal Majorana sterile 
neutrino line. A diagrammatic description of eq. (fT6l) is shown in Fig. [0 Thus, the 
smallness of this Majorana mass is natural. 

Masses m^ and Mr are determined by model parameters and then we can use 
Takagi diagonalization with the unitary matrix Vo to diagonalize the symmetric 
complex matrix [26|. The PMNS matrix Vo can be parametrized 

as 


/ 


C 12 C 13 


Vb = 


— S12C23 — C12S13S23 


V S 12 S 23 — C 12 S 13 C 23 


s] 2 ci3 ^ s; 3 ^ 

C 12 C 23 — S] 2 s 13 23 C 13 S 23 

— C 12 S 23 — s] 2 Si 3 C 23 C 13 C 23 J 


(17) 
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<*> (H) (H) (<f>) 

• • • • 

\ / \ ✓ 



Figure 1: A diagrammatic description of eq. (fl6l) 


where c tJ = cos Q i3 , .s t j = e l6 ’ :j sin 9 lJ , O^s are mixing angles of active neutrino and 
SijS are CP phase angles(3 for Majorana neutrinos). 

As we know, to suppress the neutrino masses to sub-eV with the seesaw mech¬ 
anism, we need extremely large Mrs. Interestingly, this can be achieved by the 
exponentially increasing behavior of the VEV ($(?/)). The matrix element M R lJ 
can be estimated as follows 


Mi 


R,ij 


y 




Li -1 




dy{$(y)) 2 f u m ( y)f v m (y) 

v iR U jR 

rL™ 

y {rn) N-n 2 Sij / dyv 2 cosh 2 (M(y — yo))e~ 2MN ^ v ~ 

s y {m) \Q\ \m 2 n , , m 2 n Wosr/r( ,) 

[Jr - 1 + - y o) ) 

=inh[2M(L,n - y 0 )) + ((A 

Ml 




N 


(18) 


We plot the 3rd element of the diagonal, Mr, 33 , as a function in terms of Mjv and 
let L^ } —» L, = 0.65L, 0.7L, 0.75L in Fig. [2] This function increases when 

increases or Mn decreases, and we find that if Mn < 15, ~ 0.75 then Mr , 33 can 

be as large as 500L” 1 ~ 10000L” 1 . 

Apparently there are hierarchies M^ n <C Mr , 2 2 <C Mr , 33 and one may worry 
that some element of matrix is not suppressed by Mr, 33 , but by 

Mr, n instead. So we show the explicit expression of M R 1 (m^) T as follows 


/ J^L. + dkl2_ 1 ^13_ 

Mu,11 Mu,22 Mu, 33 

mnm2i 1 1111211122 _i_ 1111311123 

Mu,11 Mu,22 Mu,33 

mnm 3 i , m 1211132 , 1111311133 

\ M R , 11 M R ,22 M R ,33 


mnm 2 i 1 1111211122 1 mi 3 m 23 

Mu,11 Mu,22 Mu,33 

m 21 1 m 22 1 m 23 

Mu,11 Mu,22 Mu,33 

1111111131 , 1111211132 1 1111371133 

Mu,11 Mu,22 Mu,33 


IIIIIHI 3 I 1 
Mu,11 “ r 
mnrrt3i , 

Mu,11 T " 
m li 1 
Mu,11 


ffl 12 ffl32 I 7771377733 
Mr ,22 Mr, 33 

7771277732 _i_ 7771377733 

Mr ,22 Mr, 33 

m 32 I m j 3 

Mu,22 Mu, 33 


\ 

/ 


(19) 
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M R , 33 K/ m) L- 1 ) 



Figure 2: M R33 vs. M N with —)■ 1 and = 0.65,0.7,0.75 


Then we see that all terms contain m 33 (which assumed to be the largest element 
of Dirac mass matrix) are suppressed by Mr i33 . Also note that mu, 77122 , etc. are 
usually much smaller than m 33 , so their suppression don’t need masses as large as 

In conclusion, thanks to the exponential-like VEV of the scalar, although our 
scale L _1 is only about order of TeV, it is still possible to lower the neutrino mass 
to sub-eV with the Majorana mass Mr. 


3.2 Numerical results and discussion 


Since we have fitted the parameters of the scalar $ and H in the quark case (see 
AppeiididDl) . we set them fixed in the lepton fitting. Although we extend the gauge 
group in this model, but it will not affect the parameters we obtained in the quark 
case. Note that the parameter y^ comes into the fitting only in a combination 
so we will not treat y (m ' ) and y ^ separately. In our fitting, we only consider 

vV m) 

the normal hierarchy of neutrino mass. 

The recent experiment data of leptons have been used in our fitting are listed in 
the following, 


• Masses of charged leptons: m e = (0.510998928 ± 1.1 x 10 8 )MeV, = 
(105.6583715 ± 3.5 x 10" 6 )MeV, m T = (1776.82 ± 0.16)MeV (Sj. 


Mass squared difference between two generations: Am 31 


10 _3 eV 2 , = (7.5 ±0.19) x lQ- 5 eV 2 23 


(2.473 ± 0.069) x 
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Mixing angles: sin 2 6* 12 = 0.302 ± 0.012, sin 2 6> 23 


0.0227 ±0.0024 23 


0.413 ± 0.032, sin 2 0 13 = 


Since there are more free parameters than data, we only show one set of the possible 
parameters. They arc listed in Table [2] If we assume that y^ ~ 0(1) (a parameter 


Table 2: Best fit parameters for leptons 


T W 

L ( l } 

T W 

M l 

0.378389 L 

0.670380 L 

0.908743L 

-11.792317 L~ x 

j-in) 


j -in) 
lj 2 

M n 

0.062289 L 

0.515437 L 

0.741436L 

13.293167 M 1 

T ( e ) 

L[ e) 

T ( e ) 

L 2 

M e 

0.317799 L 

0.448665L 

0.701578L 

36.580911 M 1 

y(e)v 

V2 

y^v / TeV 

V L 1 

- 

- 

0.317575GeV 

0.000319953GeV 

- 

- 


with a tilde means it is scaled by L to be dimensionless), then we can see that the 
hierarchy between 3^ e ' and y^ n> is about 3 order of magnitude which is acceptable. 
Notice that when L _1 has larger magnitude such as lOTeV or lOOTeV, may get 
closer to ■ If we compare the Yukawa couplings with that for the quark sector 
in Table. [4j we will find that has the same order with y( d \ So no hierarchy of 
the Yukawa couplings between quarks and leptons. All lepton 5D masses Ml,Me 
and Mtv are 0(10) up to the scale M 1 which are also seemed natural. 

This set of parameters will give 

• Masses of charged leptons: 

m e = 0.510999 MeV, m^ = 105.65837 MeV, m T = 1776.79963 MeV. 

They all deviate the experimental value less than 0.01% as the fitting required. 

• Masses of neutrinos: 

m 1 = 0.005074 eV, m 2 = 0.010092eV, m 3 = 0.049868eV. 

Comparing with the data, the mass squared differences between the 1st and 3rd 
generation deviates the experimental one about 0.5%, while the mass squared 
differences between the 1st and 2nd generation deviates the experimental one 
about 1.5%. 

• Masses of sterile neutrinos: 

M 1 = 1.2144GeV y( ^.~ 1 , M 2 = 4.9870TeVx , M 3 = 358.8498TeV^^ 

Both y ( m ) and the scale L _1 are undetermined. But we can see that if 
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i/HL -1 0(1 lOTeV), the lightest sterile neutrino can be produced by 
the LHC, and since it interact weakly with other particles, it may only con¬ 
tribute to a little part of the missing E t . 

• Mixing angles: 

sin 2 0 12 = 0.30315, sin 2 6 23 = 0.4359, sin 2 0 13 = 0.0221. 

They all deviate the experimental value less than 6%. 

• CP phases: 

S 12 = 0.1944, S 23 = 1.2796, <5i 3 = 3.0716. 


We can also calculate the effective Majorana mass as: 

(mpp) = | y, m k U'e k \ = 7,43meV 


( 20 ) 


This quantity is related to the double-beta decay which now have limit (m^) < 
120 ~ 250meV (90%CL.) j^]. Not surprisingly that our result is far from the 
experimental limit since the masses of active neutrinos are all smaller than 100 
meV. 

We can also estimate the mass of gauge held C M as follows: 


-mlC^C, = 


A/lj c 2 (4>) 2 C'‘C„ 


( 21 ) 


which implies 


m 2 = 


dyW* = 2§f m 


M 


L L 

- + -- I 

2 AM 


sinh(2M + 2 My 0 ) — sinh2 My 0 ) j 


~9l\Q\L 


-2 


n 2M(l—j/o) 


2 M 2 

which further leads to 


( 22 ) 


m r 


9c 


\/i5 ] l -' 

V2M 


e M(l-yo) _ ( 124 . g c )TeV 


--1 


(23) 


TeV. 

10 ~ 10000 TeV. Notice 


So for g c ta 0.1 ~ 1,L _1 « 1 ~ 100 TeV we have m c 
that there is another mixing effect if H is U( 1)' charged. When EW symmetry 
breaks, there will be a mass term involving Z and C [19|, then to obtain the mass 
eigenvalues we shall diagonalize a mass matrix in (Z,C) basis as 


M = 


m z 

/3m 


/3m 


mz 


(24) 
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where /3 is a factor about 0(1) or less. Since our m 2 c is apparently much larger than 


m 


Zi 


so the mixing would not be significant and the p = —f is very closed to 1, 
where m i is the smaller mass eigenvalue. Notice that this heavy gauge field will 
also significantly suppressed the effective coupling of some process mediated by it. 
The effective coupling which is similar to the Fermi constant G c ~ ^ = ^ V2iL - ij 2 ~ 
yooL-y/TcV) 2 * s muc ^ smaller than Gp, so this process will not change the whole 
amplitude. 

Interestingly, given the parameters shown in Tables. [21 [4j we do not need to 


worry about the constraints from the proton decay. Following the analysis of [15 


the dimension-eight operators lead to proton decay are QQQL , DUQL , UDEU and 
QQUE. We show the domains of the first generation wave functions which involved 
in the operators in Fig. [H]. We find that for each operator, there are at least two 
domains do not overlap, and thus the integration vanishes. 


,(0) 

h,z 


, -^-• 


d w 


* (0) 


.0 

0.2 

0.4 0.6 0.8 1.C 

y/L 

d w 



„(0) | 



U l.R 1 



* : 


J 0) 







U ,R 


fl <0) 


,( 0 ) 

l l,L 


y/L 




r~i 


1.0 0.0 0.2 0.4 0.6 

y/L 


0.8 1.0 


Figure 3: The domains of the first generation wave functions. The left-top is for the 
operator QQQL ; the right-top is for the DUQL ; the left-bottom is for the UDEU ; 
while the right-bottom is for the QQUE. 
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4 Summary 


In this paper, we have discussed the possibility to generalize the model constructed 
in Ref. 14|, [l5| to a Majorana neutrino case. The extra dimension scale L -1 is about 


several TeV, which seems far from the scale for seesaw mechanism and is unlikely 
to explain the small neutrino masses naturally. But we note that the smallness of 
neutrino masses can be a synthesized effect of the Type-I seesaw and the overlap 
integration of the localized lepton wave functions. We find that a 5D scalar $ 
with an exponentially warped VEV, which was initially introduced in Refs. 


14, 


15 


to generate a hierarchy between generations, can also be used to generate large 
Majorana masses for the neutrino right-handed 0-modes. The strategy is to let $ 
couple with singlet neutrino field in the manner ^* 2 N^Nr. When <1> acquires a 
non-zero vacuum expectation value, (Q(y)) 2 , which exponentially depends on the 
extra dimension coordinate y, will be extremely large near y = L so that the third 
generation of right-handed neutrino will be very heavy and turn on the seesaw 
mechanism. At the same time, if the positions of the O-thickness branes and the 5D 
bulk mass M^ are properly chosen, the overlap integration of the left handed and 
right handed neutrino wave functions will be also smaller than that of the charged 
leptons. Both of these effects work together, and they can significantly suppress the 
neutrino masses. 

To justify the model, it is necessary to add a U( 1)' gauge symmetry into the 
model. This symmetry prohibits some troublesome terms like Lcr 2 H* a 2 L c and the 
explicit Majorana terms. When $ obtains a non-zero vacuum expectation value, the 
U(l)' symmetry will break spontaneously. Since the mass of the U( l) 7 gauge boson 
is very large, it will not change the prediction significantly. For consistency, we also 
discuss how the anomaly cancellation conditions constrain the 17(1)' charge of each 
field. The numerical results of our model parameters have no significant hierarchy 
among them. They can fit all masses and flavor mixing data very well. We use 
this set of parameters to calculate some observable quantities such as the effective 
Majorana mass, and we find it is consistent with the double-beta decay experiments. 
Our parameters also rescue us from the stringent proton-decay constraint on the cut¬ 
off scale. 
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A The general setup of the framework 


In this appendix we briefly review the extra dimension model with point interactions. 
The basic setup is to let all fields live in 5D spacetime and have point interactions 
with some O-thickness branes HQ. The point interaction means a (5-function- 
potential-like interaction which vanishes everywhere except at a point in the 5th 

dimension BEJEL 

The action of a 5D fermion field T(x,i/) is given by 14 


S — I (ft i I dy<t(x,y)(ir M d u + M F )'b(x, y) 


(25) 


where Mp is the 5D bulk mass, and the T matrices obey the Clifford algebra 
{TjvfjTjv} = —2t}mn with the 5D metric r] MN = diag{—1,1,1,1,1} and the in¬ 
dices M, N = 0,1, 2, 3, 5 and /z, v — 0,1, 2, 3. An explicit representation of the T 
matrices is T M = and T y = —*75 = 7 ° 7 1 7 2 7 3 . The variation of the action (127>j) is: 


6S = 


<fx / dy [S^(iT M d M + M F )T + 4 <{iT M d M + M F )hT] 


d A x / dy 5^{iT M d M + M F )T - T(7T M d M — M F )<5T + 


M 


tr 


iM j 


(26) 


Thus, 5S/5ty = 0 implies the equation of motion (EOM) for T: 


(■ iV M d M + M F )T 


/'-dy + M f ia^dy \ ( T L \ _ 

V id^dy d y + M F ) 


(27) 


where the field T(a;,i/) has been decomposed into the left-handed and right-handed 
components 4 ’l,r — Pl,rS& — [(1 -F75)/2]T in the chiral representation of Dirac ma¬ 
trices 7 O Taking complex conjugate of eq. (1771) gives the EOM for T: T(iT A/ d m ~ 
M f ) = 0. Substituting it and (127]) into (1HT|) and taking 5S = 0, we obtain 


0 


d A x 


dyd M (^T M S^) 


d A x 


dy[d^T y S^) + d y (VT v 5V)] 


(28) 


Since the integral of the 4D total divergence vanishes: J d^xd^^T^d^) = 0 , we 
have 


j dyd y {^T y S^) = 0 


(29) 
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which, as we have seen, is required for the consistency of the EOMs for 'h and 
Now let us consider a toy model, in which the extra 1-dimensional space is an 
interval with length L and in the 5th dimension there are 3 boundary points assigned 
as 0, Li(< L ), L, respectively. In this case, eq. (T29|) implies 

+ [ ) dyd v (Vr v 6V) 

JLx+eJ 


0 = 


dyd v (Vr v 5V) = 


= (vr y 6v) 


- (yr y 5^) 


y=L 


(^r y 8^) 


y =o 


- (^T y 5^>) 


y=Li-e 


(30) 


y=Li+e 


where e is a positive infinitesimal length. A sufficient condition to satisfy eq. (1HU1) is 
to let the term vanish at all the boundary points: 


VT V 6V = L - 'L^r) =0 (at y = 0, L x ± 6, L) (31) 


It is sufficient to satisfy eq. (l3Tj) by imposing the Dirichlet boundary condition 


'I'r = 0 or 'I'l = 0 (at y = 0, L x ± e, L) (32) 


More specifically, we can take '&# = 0 (or 'Pl = 0) at all the boundary points to 
realize the left-handed (or right-handed) fermions in the zero mode sector, as we 
will discuss later. 

Multiplying the operator {iY N 8^ — M F ) on eq. (127|) from the left gives 


(iT N 8 N - M F )(iT M d M + M f )^> 


^—DD^ + dyd 11 


-DW + dyd 11 



where D = d y + M F , D t = — d y + M F , and 8^8^ = r] fJ , u d y, d l/ = —8% + V 2 with the 4D 
metric r) lw = diag(—1,1,1,1). Let us separate variables of the solutions of eq. (133)) 
as follows 


^L(x,y) = ^2^ ( F \x)f^ ) {y) , ^ R (x,y) = ^2^\x)f^)(y) (34) 

n n 

For every particular solution of the left-handed wave-function, F (x,y) = ^\x)f^ ) {y) ) 
we have 


0 = 


-DD' + d^)^\x)f, M (y) 


-DDtftn, (y) t/iS’ ,, (a:)+ I>l\x) Im(v) 


( n ) I 


-DD ] + M 


■^(«) i f^(y) 






X) 


where we have used the 4D Klein-Gordon equation (8^8^ — M 2 m )^\x) 
Eq. (T35|) implies 

DD^ f^(n)(y) = M^( n ) /^(») (y) 


(35) 
= 0. 
(36a) 
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Likewise, using {d^d ,J ' — = 0, we obtain 

D^Df ,(.n)(y) = M^( n ) f ,(n) (y) 

* R ' R 


(36b) 


In eqs. (I36a,p and (I36bl) . we have used the fact that the operators DD^ and D^D 
are supersymmetric quantum mechanical partners HQ I 22 I and thus they have 
exactly the same eigenvalues except for the lowest zero eigenvalue. It can be easily 
explained as follows. If f^ n )(y) is the eigenfunction of DD^ with the eigenvalue 
M2 ( „) and Mj (n) ^ 0, then 


D*D 


>(v) 


= £>> 


DD'f M (y) 


= MJ W BV wfe) 


that is, f ^ n ) (y) is an eigenfunction of D^D with the same eigenvalue M| (n) . 

f^n){y) oc D*f M (y) and let f^(y) have the same normalization as /^(y): 


(37) 


Dehne 


(Lw(y)\Lw(y)) = dy\f (n )(y)] / (n)(y) = l 


which implies 


dy \D ] f ,(n)(y)l D^f ( n ) (y) = [ dy [/.(„)(y)l DD*f M (y) = Af; 


/ / \ ,L \R ) I 1 v • “) \ rj l I / / V 

'Pi 'J j » i/n 

Then it is sufficient to get (/m»)(y)|/m«)(y)) = 1 by letting 


'0( n ) 


/^(l l )(y) — 71 r (y) 


Af.,.^ ^ 


Multiplying the operator D on the above equation from the left gives 

/^(»> (y) = — Df qh (n)(y) 


M.l. (n) ^R 


rL ' 

Substituting a pair of chiral modes of (134|) into eq. (127)1 . 

1X 71 ) 1 


we have 


/ £)t zcffidA M (®)/^)(j/)' = 
J \^\x)f^ n) (y) 


l \ x ) [^/^(y)! + Q^^Oe)] ./>)(?/) =0 


ia^d^ix) f M (y) +^‘\x) Df (n) (y) 


( n )i 


= 0 


(38) 


(39) 


(40a) 


(40b) 


(41) 


(42a) 

(42b) 
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which, together with eqs. (I40ap and (I40bj) . lead to 


(43a) 

(43b) 


ia^d^\x) + = 0 

ia^d^^ix) + +M^(n)'ifjR ) (x) = 0 


that is, 

(M^n) icr^dA ftp l\ x )\ ={] 
Wdp \y%\x))- 


(44) 


Thus, the combination ^ n \x) = ytp^\x), iPr\x)j obeys the 4D Dirac equation 

(i$ + M^(n) ) '0 1 ' n \x ) = 0 and forms a Dirac spinor. 

Suppose that the eigenequation (I36al) of DD^ has a zero eigenvalue M^ (n) = 0 
with the corresponding eigenfunction / ,(o) ( y ) called the 0-mode. That is, DD^ f (o) ( y ) 

VL TL 

0. It is sufficient to satisfy the above relation if f Ao)(y) is annihilated by D': 


= (-a„ + M F )f^(y) = 0 (45) 

• If the Dirichlet boundary condition T i (a;, t/) = 0 is imposed at y — 0, Lxie, L, 
that is, 

4(o) (?/) = 0 (at y = 0, Li ± e, L) (46) 

then eqs. (|45l) and (1461) imply that f xo)(y) = 0 at all points. Thus, the 0-mode 
eigenfunction of DD^ does not exist in the boundary condition of (1461) . 


• If the Dirichlet boundary condition T#(x, y) = 0 is imposed at y — 0, Li±e, L, 
that is, 

4<o)(j/) = 0 (at y = 0, Li ± e, L) (47) 

then this boundary condition has no effect on the equation (145|) , but the setup 
of the 0-thickness branes’ positions itself can split the solutons of (1461) into 
two independent degenerate modes: 


f N ie M ^ 
,<■)(») = \ o 

hl 0) ,(2)^ — { N 2 e Mpy 


(o <y< Lx) 
(Lx < y < L ) 

(0 <y< Lx) 
(Lx<y < L) 


(48a) 

(48b) 


where Nx and N 2 are normalization constants and, by using (1381) . they can be 
figured out as 


N 1 = 


2M7 


^2MpLi _ ^ 5 


No = e - MpLl 


2 Me 




(49) 
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Using the Heaviside step function 9(y), we can also write the two degenerate 
zero modes as follows 


t/,™,(1)(») = \f^Sf-j eMn -»)] < 50a ) 

4<»\ -i eMF ^ -»)i ( 50b ) 

The 5D wavefunction of 0-mode may be expanded with respect to 

and )(y ) as 

^ l\ x , y) = ^il( x ) (!)(?/) + ( 5 i ) 

where the coefficients (x) and i/j^l ( x ) are identihed with the 4D wavefunc- 
tions of two generations of left-handed fermions in this toy model. 

Likewise, consider the 0-mode eigenfunction f ao)(v) of D^D. It obeys the equation 
D^Df (o) (j/) = 0. A sufficient condition of this equation is 

Df(») = (d, + M F )f 4 „> (y) = 0 (52) 

• If the Dirichlet boundary condition (1461) for the left-handed fermion is imposed, 
then it is the location of the point-interaction positions, rather than eq. (I46T) . 
that affects the solutions of ()52l) and splits them into two degenerate modes: 

4™,(l)fe) = - »)] (53a) 

/<>, [% - L,)0(L - y)] (53b) 

The expansion of the 5D wavefunction of 0-mode y) with respect to 

the two modes is given by 

^r\ x ,v) = ^iR( x )f^°\ { i)(y) + i>2R( x )f^<g\(2)(y) ( 54 ) 

where the 4D wavefunctions 'ipin (x) and ( x ) belong to two generations of 
right-handed fermions in this toy model. 

• If the Dirichlet boundary condition (1471) for the right-handed fermion is im¬ 
posed, then eqs. (I52|i and (1471) imply that f,(o)(y) = 0 at all points. That is, 
the 0-mode eigenfunction of D^D vanishes in this boundary condition. 
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To sum up, if the boundary condition L = 0 is imposed at all the O-thickness 
branes’ positions, then the 5D fermion held ty(x,y) has only right-handed 0-modes 
(x, y) as given in eq. (I5T|) : instead, if R = 0 is imposed at all the boundary 
points, then ^(x,y) has only left-handed 0-modes ^^\x,y) as given in eq. flSTJ) . In 
a word, the Dirichlet boundary condition ^ l,r — 0 makes the 0-mode wavefunctions 
of ^(x,y) to be chiral. Including the KK modes (he. the modes with M^ (n) ^ 0), 
the expansion of a 5D fermion held T(x,|/) in all modes is given by 

• For \lR — 0 at y — 0, L 1; L 

v( x ,y) = \j x _ 2 ^2M F L; e ~ MFy \. e (y) e ( L ^ - y)]^lfl(^) 

+ \ h _ e~ MHy ~ Ll) [8(y ~ U)0{L - y)]^(x) 

+ (KK modes) (55) 


For = 0 at y = 0, L 1; L 
^( x ,y) = 


2 Mp 


^2MpL\ _ ^ 


e MFV [e(y)6(L 1 - y)}^{x) 


+ 


2 Mp 


(KK modes) 


e"' 1 ' Ll) [S(y - L{)e(L - y)]i!>fl(x 


(56) 


To realize both left-handed and right-handed 0-mode fermions in this 2-generation 
toy model, we need at least two 5D fermion holds, ^i(x,y) and ty 2 ( x ,y)- O ne 
5D fermion T 1 (a:, y) has two left-handed 0-modes due to the boundary condition 
PrT i(x,y) = 0 at points y — 0, L i, L; while another 5D fermion ty 2 (x,y) has 
two right-handed 0-modes from the boundary condition P F ^ 2 (x,y) = 0 at points 
y = 0, L[, L. The locations of and L[ are in general not equal. Indeed, it is the 
inequality of L\ and L\ that leads to the mixing of the two generations of fermions. 
A schematic picture of the wave functions of these 0-mode chiral fermions is shown in 
Fig. [4j To give the chiral fermions masses, we need to introduce an extra 5D scalar 
held <F(x, y), which will acquire a nonzero VEV after the electroweak symmetry 
breaking. The mixing structure of the Dirac mass matrix is also explained in Fig. SJ 
In addition, it is worthy to point out that the operators D = d y + M F and 
D t = -dy + M f can be used to construct a pair of supersymmetric generators, 
Q = D^°Pl and Q I = D'^Pr, which satisfy the supersymmetric algebra (See the 


paragraphs between eqs. (8) and (9) in Ref. [llj for more details): 


Q 2 = Q ]2 = 0, {Q,Q*} = 2H, [Q,H] = [Q\H] = 0 


(57) 
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Figure 4: A schematic diagram of wave functions for chiral 0-mode fermions. The 
red curves represent the wave functions for two generations of right-handed 0-mode 
fermions, while the blue curves represent the wave functions for two generations of 
left-handed 0-mode fermions. The black line is a profile of a scalar $’s VEV. The 
overlap integration of the profiles in different interval gives the corresponding mass 
matrix element. 


24 











The Hamiltonian operators (up to a constant factor) is H oc {Q,Q^} = DD^P R + 
D^DP l , and the pair of modes (/^(«) (y)'i^ < l\x), f^(n)(y)i^ R \x)) T is an eigenstate 
of H with eigenvalue M^ (n) . 


B Quark masses hierarchy and flavor mixings 


The Yukawa terms which generate the masses for quarks are: 

Cltks = ~J dy{yM$Q(ia 2 H*)U R + yV>VLHD R + h.c] (58) 


where y^ and y^ are the couplings with dimension —2 for the up type and down 
type quarks, respectively. 

Note that we will let $ to be £/(1)' charged. Then if we don’t want the t/(l)' 
breaks explicitly, we should also make U R , D R , Q and H to be U( 1)' charged. We 
have determined the t/(l)' charge for each field in section [3] We can see that terms 
as Q(ia 2 H*)U R and LHD R can be forbidden by the U{1)' symmetry. 

After the U (1)' and electro-weak symmetry breaking,we obtain Dirac mass terms 
of quarks. The mixing structure of the mass matrix will be generated by the overlaps 
of wave functions from different generations. Then we can write down the mass 
matrices as 


m^ u> = 


/ 

m 12 

777^^ 

| , m «> = | 

/ mfj 

mf 2 

™1 3 \ 

° 

777-2 

m 23 

0 

777^ 

77723 

V 0 

0 

rri-k) 


V o 

0 

miJ 




mif =y (d) I dyf a m(y)fj»(ymy))(H{y)) 


dyfm (y)f u m (y) (<j>(y)) (H(y)*) 

ViL a jR 




“j'H 


(59) 

(60) 
(61) 


The integration range (a, b ) represents the overlap region between the profiles / (o) (y) 

QiL 

and / (o )(y) or f,(o)(y). The integration will contribute to a diagonal element when 

u jL a jL 

i = j, and an off diagonal element when i ^ j. Two Dirac mass matrices and 
m (rf) are apparently complex and we can diagonalize them with unitary matrices 
h“>(h i) ) and v£>( V^). 


AS, = 


(62) 


Then we can compare the masses with experimental data. Using matrices V\' a) and 
V[ d \ we can calculate the CKM matrix which is defined as 

Vckm = (63) 
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The CKM matrix contains not only information about flavor mixing angles but also 
information about the CP violation. The CP violation can be characterized by the 
Jarlskog invariant J defined as 

3 

bn [{ycKM)ij{VcKM)kl{V^ KM )u{V^ KM )kj\ = J ^ tikmtjln (64) 

m,n =1 


We list the experimental data used in our fitting as follows 

• The up and down type quark masses are shown in Table. [3] 


Table 3: Quark masses from ref. 25 


up type quark 

mass 

down type quark 

mass 

u 

2.3 ± 0.6 MeV 

d 

4.8 ± 0.5 MeV 

c 

1.275 ± 0.025 GeV 

s 

95 ± 5 MeV 

t 

173.5 ± 1.4 GeV 

b 

4.18 ±0.03 GeV 


• The absolute values of CKM matrix elements from ref. 



are 


| VcKM 


/0.97425 ± 0.00022 
0.230 ±0.011 
\ 0.0084 ±0.0006 


0.2252 ± 0.0009 
1.006 ±0.023 
0.0429 ± 0.0026 


0.00415 ± 0.049X 
0.0409 ±0.0011 
0.89 ±0.07 ) 


(65) 


The Jarlskog invariant from ref. 25 


is J = (2.96 ±0.18) x 10" 5 . 


After fitting the data listed above. We find a set of parameters, which is compatible 
with the data, and show them in Table. SI We have set |A[ = \XL\ = 0.001, \Q\ = 
\QL 5 \ = 0.001 and y 0 = y 0 L~ l = —0.16 fixed as reference |l4| did, so the only free 
parameter of $ is M. Since $ and H also couple to leptons,the values of M and 
6 which are found in the quark case will be set fixed to reduce the number of free 
parameters in the lepton case. In the following,a parameter with a tilde means it 
has been scaled to dimensionless by multiply some power of L. 

Note that we can calculate L± in the Robin boundary condition by: 


L + = = —0-118074 L 

L - = M)= °-l° 4502i 


( 66 ) 


Then we find that M = 9.36099 < j- — 9.5692, which is consistent with the 
symmetry breaking condition \M\ 2 < -pr —. 

Using the parameters of $ we can calculate the tree level mass of the 4D exci¬ 
tation 4>(x). One of its degree of freedom will be gauged out by the gauge boson 
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Table 4: Best fit parameters for quarks 


r(«) 

^0 

L? } 

r (?) 
lj 2 

Mq 

0 

0.31423L 

0.67665L 

9.26018L" 1 

T (“) 

V) 


r O) 

M v 

0.05218L 

0.06095L 

0.56328L 

—4.48152L” 1 

r O 

L 0 

L {d) 

T A) 

^2 

m d 

0.11866L 

0.23128L 

0.66636L 

5.71010 L” 1 

M 

y (u) v/V 2 

y {d) v/V 2 

6 

9.36099L” 1 

3.15684 GeV 

0.20552 GeV 

2.91684 


of U(l)' when the symmetry breaking occurs. To obtain the mass of (j>(x), we shall 
consider its excitation around the minimum of potential 

£[$]=[ dy{-&dfo + M 2 |$| 2 +^|<L| 4 } (67) 

Jo 2 

Substitute the zero mode (y)[y + 0), vf^°\y) = (<£>(?/)) into £[$] and use 

the minimized condition:— dyf 0 (y) + M 2 / 0 + Az/ 2 /q = 0, we can get the mass 

ml = 2 £ dy(2X[<Hv))yi) « = ^«(!-«„, L - 2 (68) 

which implies 

m 0 « V^i e M(i-yo) L -i ~ 5.55TeV • ^ (69) 


If the scale L" 1 ~0(lTeV), this mass is under the energy scale of LHC. But it is 
unlikely to be detected in the recent experiments, because the the 0-fermion-fermion 
couplings are so weak. This can be seen by estimate the couplings as 


where 


Aq) 
’ij 


m. 


(<?) 

ij 


■A 


( e ) 

Ae) = m ij 


A 


v 


(70) 


A 


_ \/VI _ 

\J 1 + sinh(M) cosh(M — 2 My 0 )/M 


[2 2VM 

V L e M(i-yo) 


v ~ 


- L 2 

M 


(71) 


(72) 
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Using the parameters in our fitting, we find the Yukawa couplings for 0-quark-quark 
and 0-lepton-lepton are 


4 9) - 0.03 x 


in, 


(?) 

ij 


-1 


r e ) m\f 

Q? ^ 0.03 x 13 


-l 


(73) 


Both Yukawa couplings are much weaker than the Yukawa couplings for Higgs- 
quark-quark and Higgs-lepton-lepton. Since the coupling is proportional to the 
mass, the strongest Yukawa coupling may be the coupling of 0-top-top which is 
about 0.03 x 0.17 ~ 0.005 when L -1 ~ 1 TeV. 

Note that there is a Cj < f>| 2 |//| 2 term may lead to some problem with the gauge 
universality as discussed in reference 14J. We will just let C to be small enough 


(about 10 'for L 1 ~ 1 TeV) to resolve this. 


C Why an explicit Majorana mass term does not 
work 


The 5D charge conjugation operator C is defined as 

CT M U _1 = (T m ) t 

with properties: 

c T = c- 1 = C f = -C 


(74) 


(75) 


It is easy to check that C can be written as C = 7 ° 7 2 (* 75 ) jUJ]. We can write it in 
Weyl basis 


C = 


Cab 


— 6 


ab 


(76) 


The charge conjugation of a 5D fermion is defined as 

T c = C^> T 


We can also write it down in Weyl basis: 


Xa(x, y) 

-P°(.x,y) 


(77) 


(78) 


Note that the relation (T c ) c = T no longer holds in 5D case and the correct relation 
is (T c ) c = —T. 






Now we consider to add terms as ^>iY M 8 m^ c + h.c after several lines of calcu¬ 
lation, we can get 


WT M d M V c = ^iT M d M C^ T 


'M'~ 

^M,T r C\ 




(79) 


This implies that these terms can be absorbed into the boundary terms and do not 
contribute to the equations of motion. 

ffowever, the mass terms as Mr\E'T c + h.c. survive and will contribute to the 
equations of motion. Now let’s add the mass terms into the action: 


S 


d 4 x J dy[^f(x,y)(iT M d M + M F )V(x,y) + 1 -{M R ^ C + h.c.)] (80) 


The variation of the action flSDJ) is: 


6S 




dy 


dy 


5^>{iT M d M + + V(ir M d M + M f )5 T 

+ ^M r 6VV c + ^M R ^S^ C + ^M r SW^ + i M r W5 T 
6^(iT M d M + M f )T - T( iT M< d M ~ M F )5 T 


+ <9 M (^T M hT) + 


(81) 


Thus, the equation of motion (EOM) becomes: 

0 = {iY M d M + M F )y + M r ^ c 

f-d y + M F ia^dy \ t £ a (x,y)\ + f M R \ ( Xa(x,y) \ (82) 

V iv 1 * 9 !* dy + M F ) \x' & (x,y)J \ M r) \~£ h { x ,v)) 

If we try to separate the field in modes as £ a {x,y) = f^ n Ky)^ l \x), Xa = 
9 ^ n \y)Xa\x) , then the equations for each mode become 

(~d y + M F )fW(y)&\x) + M R g^\y) X { :\x) + g^*(y)ia^x^\x) = 0 (83) 

(d v + M F )gW*(y) X ™ & (x) - M R f^*(y)^(x) + f n) i^d^\x) = 0 (84) 

Apparently, in a general case, it is impossible to factor out the functions f^ n \y), g^ n \y) 
from the 4D Dirac equations of spinors £ a ( x ), Xa( x )■ This means A special choice 
which can achieve this is to let M F = 0 and x^ a = — £ , then the EOM become: 

(dy + M R )£a(x, y) + iahd^ a (x, y) = 0 
(~d y + M R )^ a (x, y) + id^ aa d^ a (x, y) = 0 
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(85) 

( 86 ) 






We can recover the 4D Dirac equation for a Majorana fermion by setting £ a (x, y) = 
A£ a (x) where A is a constant so the profile is independent of the 5th dimension 
coordinate y. Thus, this fermion has only one mode with a Majorana mass Mr. 
But this solution requires some special choice of the 5D fermion. 

If we accept this special pattern of fermion to be the singlet neutrino Nr, and 
generate Dirac masses with the Yukawa interaction, then the seesaw turns on when 
the Majorana mass is much larger than the Dirac ones. However, an operator as 
Ta 2 H*Wa 2 L c is still allowed and it will contribute to the Majorana masses of left- 
handed zero modes. Now we have to diagonalize the following mass matrix 


M 


Ml M d 

m t d m r 


(87) 


In the large Mr limit, the light neutrino masses are m v ~ Ml — MrMr/Mr. These 
masses should be as small as 0(0.1 eV) to fit the current neutrino mass bound and 
imply that either we use an unnaturally small coupling for La 2 H*a 2 L c operator 
or we fine-tuned the parameters to cancel Ml by MrM^/Mr in high precision. 

Actually, in the SM the gauge symmetries and the lepton number conservation 
do not allow the explicit Majorana mass term and La 2 H*H^a 2 L c to exist. However, 
in this model we are going to add a SM-gauge-group singlet neutrino field into the 
model, and try to violate the lepton number explicitly. Thus, we have to face these 
annoying terms unless they are also forbidden by some symmetry. The strategy we 
use in the paper is to forbid both La 2 H*Wa 2 L c and MR^>fy c + h.c. terms by a U(l)’ 
symmetry. Then the singlet neutrinos have chiral zero-modes as any other fermions. 
Their right-handed Majorana masses are generated by the VEV of $ with the same 
mechanism as their Dirac masses generated by the VEV of $ and the Higgs field. 
Now the mass matrix we need to diagonalize is 


M 


( 0 M d \ 

\M T D m r ) 


( 88 ) 


In large Mr limit, the light neutrinos masses are m u ~ MrMr/Mr which can be 
naturally suppressed to 0(0.1 eV). 
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